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, We evaluate analytically some ground state properties of two-dimensional harmonically confined 

■ Fermi vapors with isotropy and for an arbitrary number of closed shells. We first derive a difi'erential 

pL^ ' form of the virial theorem and an expression for the kinetic energy density in terms of the fermion 

particle density and its low-order derivatives. These results allow an explicit difi'erential equation 
to be obtained for the particle density. The equation is third-order, linear and homogeneous. We 
also obtain a relation between the turning points of kinetic energy and particle densities, and an 
expression of the non-local kinetic energy density functional. 
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! I. INTRODUCTION 

C/5 _ 

We have previously studied in detail the ground-state properties of harmonically confined non- interacting fermions 
^ ' in one dimension because of current interest in experiments on ultracold gases of the '^''K and ^Li fermionic 

isotopes populating hyperfine states inside magnetic traps The s-wave collisions between fermions in the same 

hyperfine state are suppressed by the Pauli principle, while p-wave scattering and dipole-dipole magnetic interactions 
are very weak at very low temperature jH). A one-component gas of Fermi atoms in a fully spin-polarized state 
Q ■ inside a magnetic trap thus is a close laboratory realization of the ideal non-interacting Fermi gas in an external 
O harmonic potential, allowing access to the kinetic energy functional which is invoked in the density functional theory 
I of inhomogeneous fluids |^ . 
I ' The interest in such experiments has motivated us to study by essentially analytical methods the extension of the 
^ , results of Ref. [Q and j2j to treat independent fermions with two-dimensional (2D) isotropic harmonic confinement, for 
^sD ' arbitrary number of closed shells. Throughout we study singly filled levels because of the experimental situation 

OO , involving spin-polarized fermions. 

' The outline of the paper is then as follows. In Sec. |l^ below, we derive an explic it fo rm of the kinetic energy density 
\ for 2D harmonic confinement for an arbitrary number of closed shells. In Sec. Ill a differential form of the usual 



I 



o 



integral virial theorem, but now specific to 2D harmonic confinement, is invoked and this allows the kinetic energy 
density result of Sec. |l| to be rewritten as a differential equation for the particle density. Sec. IV gives the derivation 



of a relation between the turning points of the kinetic energy and particle densities. This allows the construction 
' of an explicit form of the kinetic energy density functional solely in terms of the particle density and its low-order 
derivatives. Illustrative numerical results are presented in Sec. ^ The final section constitutes a summary, together 
with some proposals for further studies. 



II. DERIVATION OF THE RELATION BETWEEN AVERAGE KINETIC ENERGY DENSITY AND 

PARTICLE DENSITY 



, Three forms of the kinetic energy density are used in the paper: (i) t[r) = — (fi,^/2m) i/iiV^-i/ii, (ii) ta{r) = 
■ (S^/2m.)X]i I^V'iP aiid ft^'^) an average of t{r) and tcir) called t{r). The quantity 2i{r) equals the trace of the 
5^ ] tensorial momentum flux density and the difference tcir) — t{r) is simply [ti^ / Am)V'^ p{r) where p{r) is the particle 
density. All three kinetic energy densities defined above lead to the correct total kinetic energy T = / cPrt{r) etc. A 
brief reminder of the solution of the Schrodinger equation for the problem of present interest is given in Appendix 

Our starting point is the result for the 2D fermion particle density in the presence of the harmonic confinement 
V{r) = mw^r^/2, as obtained by Brack and van Zyl 0, namely 

1 ^' 

Pi^) = E(*^ + 1 - M)(-l)^ip(2/) exp(-2//2) (1) 
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after adaptation to the case of single occupancy of the trap levels. Here aho = {h/mcoY^'^, y = 2{r/aho)'^, [M + 1) is 
the number of filled shells and L^(y) are the Laguerre polynomials. The average kinetic energy density t{r) is given 

by 



M 



27ra2 — 

ho ^=0 



^(M + 1 - ,if{-lYL^{y) exp(~y/2). 



(2) 



We now define the following summations over the Laguerre polynomials: 

M 

Qi''^'\y) = E f^H~^rL,{y) : z = 0, 1, 2. 

Then the particle density and the average kinetic energy density are readily expressed in the form 

1 



p{r) 



"ho 



(M + l)Qi''+'\y) - Q^^'+'\y) exp(-y/2) 



and 



t{r) 



huj 

2^ 



(Af + l)^Qi''+'\y) ~ 2(M + \)QT^'\y) + QT^'\y) exp(-y/2) 



(M+l) , 



.(Af+1), 



The next step is to utilize the differential equation for the Laguerre polynomials, namely 



dy 



dy 



This is readily shown to lead to the following relations for the summations Q[''^^^\y)'- 



d^ , , d 



dy'^ 



dy 



Q[''^'\y) + Qil\^'\y)-o 



(3) 



(4) 



(5) 



(6) 



(7) 



for i — and 1. These allow us to express p{r) and t{r) solely in terms of Qo^^^^\y) and of its first two derivatives, 
leading after some manipulations to the result 



p{r) 



(8) 



Equation (||) is one of the central relations for the present study. 



III. DIFFERENTIAL EQUATION FOR FERMION PARTICLE DENSITY 



Having related the kinetic energy density i{r) to the particle density p{r) and the potential energy in 
Eq. (j^), for an arbitrary number of closed shells, we shall next invoke a differential form of the usual integral virial 
theorem. This we write in the form: 



dt{r) , ^dV[r) 



(9) 



For the specific case of 2D harmonic confinement with filled shells, the differential virial theorem (^) can be proven 
starting from Eqs. and (|^) for any given number of shells. The demonstration is given in Appendix |^. 
Forming the virial r • F of the force F, as —rdV{r)/dr, and integrating over 2'Krdr we find 



/ 2Trr^^^dr = (r • F) = / d'^rpMr 
Jo or J 



dV{r) 
dr 



(10) 



Integrating the LHS of Eq. (|T^) by parts we obtain 2T for the integral value, T being the total kinetic energy, which 
confirms that Eq. (||) contains the integral virial theorem. 
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Having established Eq. (H), we can combine it with Eq. relating t(r) to V{r) and p(r) to obtain a a third-order, 
linear, homogeneous differential equation for the particle density p{r) for (M + 1) closed shells: 



8m dr 



[VV(r)] 



M 



huj 



dp{r) 
dr 



+ muP'rp{r) — 



As one immediate check of Eq. (0), it is readily verified that 

p(r) =AAexp(-rVaL) 



(11) 



(12) 



where A/" is a known (irrelevant for present purposes) normalization factor, satisfies Eq. ( pi] ) for M = 0. 

Equation ( pi] ) achieves a long-term aim of density functional theory for the present case of 2D harmonic con- 
finement of independent fermions: namely to allow the direct calculation of the particle density from the potential 
energy {muj'^r'^ /2 in this case), without recourse to individual wave functions. It constitutes the two-dimensional 
generalization of the result of Lawes and March |^ for one-dimensional (ID) harmonic confinement of fermions. 



IV. KINETIC ENERGY FUNCTIONAL 



The first object of the present section is to obtain a relation between t'ir) = dt{r)/dr and p'{r) = dp{r)/dr for this 
isotropic 2D harmonically confined fermion assembly. The result for ID harmonic confinement is (see e.g. 0) 



^ Nhuj 



We first use the definition of t{r) and the relation to{r) = t{r) + (?i^/4m)V^p(r) to obtain 

ft2 



t(r) = t{r) 



8m 



y'pir) 



(13) 



(14) 



After differentiation of Eq. (|lj) with respect to r, the use of the virial theorem 
for p{r) allows us to obtain 



and of the differential equation (hj 



p'{r) 



M - 



9 2 

muj r 



(15) 



which is the desired generalization of Eq. (fi3| ) to 2D harmonic confinement. Evidently Eq. (p5| ) relates directly the 
turning points of kinetic energy and particle densities for any chosen number of closed shells. 

Equation ( p5| ) allows us to eliminate the presence of the external potential and to obtain the non-local kinetic 
energy density functional solely in terms of the particle density and of its low-order derivatives. To this aim, we first 
write an equation for the quantity Q{r) = t'{r)l p'{r) by differentiation of Eq. ( p^ ) and by using Eqs. (|) and (||): 



d_ 

dr 



Qir) 



p{r) 



1 d 



8m p'^{r) dr 



V'p{r) 



As a second step, we obtain the function Q{r) by integration 



Qir) 



P(0) 



P(^) 



(16) 



(17) 



the integration constant being fixed by the condition Q(0) = (M -I- 3/2)huj. A further integration and elimination of 
the high-order derivatives through integration by parts leads to the final expression 



t(r) 



twir) 



16m Jo P'^(s) 



V(f) 
P(s) 



(18) 



where twir) = {fi^ /8m)[p' (r)]'^ / p{r) is the von Weizsacker "surface" contribution to the kinetic energy density j|] and 
the constant C is given by C = (M + 3/2)fia;/{2p{0)) — {fi^ /16m){'V^ p{r))r=o/ p^ iO) . This provides a generahzation of 
the result obtained by March et al. to the 2D harmonic confinement. By integration of Eq. (^8|) over the system 
one obtains an explicit expression for the kinetic energy functional. 
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V. ILLUSTRATIVE NUMERICAL RESULTS 



Equations and (g) can be employed to numerically evaluate all the ground state properties defined in Sec. ||. 

In Fig. |l| we plot the exact kinetic energy densities t{r) and tc{r) for 20 filled shells (M=19). There is seen to be 
a tiny negative region of t{r), whereas tQ{r) is by definition everywhere positive. In the same figure we also show the 
von Weizsacker kinetic energy density tw{r). It is seen that twir) is practically zero until the classically forbidden 
region is reached, while well outside the classically allowed region it eventually merges with tdr): i.e. tQ{r) — > tw{r) 
for r greater than the Fermi radius. This property is shown in the inset of Fig. |l|. 

In Fig. 1^ the differential virial theorem is illustrated in the case of harmonic confinement by comparing the 
negative of di/dr with the quantity •muP'rp{r) for various numbers of filled shells. 



VI. SUMMARY AND PROPOSED FUTURE DIRECTIONS 



We have derived a set of analytic relations for the ground state properties of a 2D Fermi gas under harmonic 
confinement, with the proviso that the Fermi particles always fill (M + 1) closed shells. Our work was aimed mainly 
at transforming the description of the system from one given in terms of wave functions to one which takes as the 
main variable the particle density. This procedure, typical of Density Functional Theory (DFT), is here possible in 
an explicit way due to the simplicity of the system, but may serve as a guideline to the use of DFT in more complex 
systems. 

Five basic equations have been obtained, namely Eqs. 



M): (|15D and (|18|). Equation (B) is a differential 

relation which gives the average kinetic energy density in terms of the particle density, and Eq. ( D ) is a differential 
version of the usual virial theorem. These relations allow the derivation of Eq. (11), a third-order differential equation 
for the particle density profile as determined by the external confining potential and by the number of filled shells, 
and Eq. ( [1^ ) , a relation between the turning points of the kinetic energy density and the particle density. The latter 
is finally used to obtain the kinetic energy density functional ( [l^ ) in terms of the particle density, which is a main 
task in DFT. 

What seems to us therefore to remain outstanding in 2D is to solve the differential equation ([Tl]), hopefully ana- 
lytically and without summation over shells. The motivation for expecting this is the "folklore" that the properties 
of a fermion ffuid are determined by its highest occupied level ||]. Of course, the challenge of the three-dimensional 
oscillator problem remains as a high priority for further study. 
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APPENDIX A: TWO-DIMENSIONAL FERMI GAS IN HARMONIC CONFINEMENT: WAVE 

FUNCTIONS AND DIRAC DENSITY MATRIX 



The Schrodinger equation in 2D reads 

. . 1 



— V^V + -mu?r'^^lj = e^} . (Al) 
2m 2 



Using afio = (ft/mw)^/^ and Huj as units of length and energy, the eigenfunctions of Eq. (Al) are expressed through 
the Laguerre polynomials L5"(r^) as 

^n,mir, <P) = , • r"^K{r^) exp(-r2) exp(zm<^) , (A2) 

\i:(n -\- my. J 
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in correspondence to the energy eigenvalues 



= 2n + m + 1 



(A3) 



The various shells are filled by taking all integer values of n and m which satisfy 2n + m = A^, N being the shell 
index. 

The "off-diagonal density", that is the Dirac density matrix can be immediately written from Eq. (A2) as 

M N 



n(r, r') 



1 



■ exp 



iV=0 h=Q 



[N -h)\ 

h\ 



(rr')2'-^i^^-f (r'^) exp[z(2/. - N){<t> - ^')] ■ (A4) 



From this matrix, the density profile p{r) for [M + 1) filled shells is obtained by taking r' — r, the result being the 
same as the one obtained by Brack and van Zyl Similarly, the kinetic energy density can be directly constructed 
in either the form of t{r) or tQ{r) by appropriate differentiation on the Dirac matrix (A4), and afterwards allowing r' 
to tend to r. 



APPENDIX B: PROOF OF THE DIFFERENTIAL VIRIAL THEOREM FOR A 2D HARMONICALLY 

CONFINED FERMI GAS 



We demonstrate here Eq. (||) for the specific case V{r) = ma'^r^/2, which is the one relevant for the present work. 
The units ato and hw for length and energy are used throughout. 

As a first step, we derive from Eqs. (]l|) and Eq. (||) recursion relations for the particle density pM+i(r) and the 
average kinetic energy density tM+i{r) of a system with [M + 1) filled shells in terms of those of a system with M 
filled shells. These relations are 



and 



PM+i{r) = PM{r) + Q^'''+^)(y)exp(-y/2)/7r 



tM+i{r) = hiir) + pM{r) + Q\^'^^\v) exp(-y/2)/(2^) 



(Bl) 



(B2) 



These involve the function Q\^^^^\y) defined in Eq. (||). 

We then proceed by induction. The differential virial theorem is readily proven for the case M = 1 of a singly 
occupied shell by direct substitution of the expressions pi{r) = exp(— r^)/7r and ii{r) = exp(— r^)/27r into Eq. (|S|). 
Assuming that the theorem holds for M filled shells, we are left to prove that 



d_ 

dr 



PM+i{r) - 



27r 



TT 



(B3) 



Physically this means that the theorem holds separately for each filled shell. 

By the use of Eq. (|) for the density profile, Eq. ^ can be rewritten solely in terms of the functions Q[''^^^^\y) 
with i = 0, 1: 



dy dy 
We can now employ the definition of the Laguerre polynomials 



k=0 



i-yy 

k\ 



(B4) 



(B5) 



and of the functions Qf'^^^\y) to transform Eq. (B4) into a double sum: 



M 

5^(-ir5^(2M+l-2M) 



fc=0 



ki 



M 



-^(-ir^(M-M) 

fj,=0 k=0 



^J■ 

p — k 



k\ 



(B6) 
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Exchanging the order of the sums allows us to collect all terms of the same order in k: each coefficient must then 
vanish independently because of the principle of identity between polynomials. Thus, we have finally to prove the 
simplified expression 



M-l 



i-ir [{2M 1 2/i)(M + 1) - (M - /i)(fc + 1)] = 



(B7) 



The above finite sum must be shown to yield zero for any integer value M and k < [M — f ). Thi s statement is proven 
by applying twice again the reasoning by induction on the index M: the validity of Eq. (B7) for M = 1 is easily 
checked, and assuming that it holds for a fixed M leads to prove that 



A/-1 



Y^{-\r{2^Ji~k+l] 



fi—k 



= (-1) 



M+l 



M! 



{M -k- 1) 



(B8) 



Applying again induction to the LHS of Eq. ( B8 ) leads to the proof of its validity, which implies the validity of the 
differential virial theorem. 
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FIG. 1. Kinetic energy densities (in units hijj/a\^) as given by the different definitions in Sees. ^ and ^ for a 2D Fermi gas 
in isotropic fiarmonic confinement witli 20 filled shells, as functions of the radial coordinate r/aho- Solid line; tair), dotted 
line: t{r), dashed line: tw{r)- The inset shows an enlargement of the turning point region, in the same units. 




FIG. 2. Illustration of the differential virial theorem: negative of di{r)/dr (solid line) and Tnuj^rp{r) (dots) in units Tuo/af^^ 
as functions of the radial coordinate r/aho, for a 2D Fermi gas in isotropic harmonic confinement with 20, 15 and 10 filled shells 
(from top to bottom). 
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